
Subroutine InitEphemTimeRangeDE()
    implicit none
    REAL*8 EthenTimeRange(3)
    COMMON /EthenTimeRange/EthenTimeRange

    CHARACTER*6  NAMS(400)
    DOUBLE PRECISION  VALS(400)
    INTEGER  NVS
    !--------------------------------
    CALL  Const (NAMS, VALS, EthenTimeRange, NVS)
end

! =========================================================

subroutine InitGMsDE()
    implicit none

    real*8 GM(11)
    COMMON /GM/GM

    real*8 kmAE
    common /kmAE/kmAE

    CHARACTER*6  NeedNames(11)
    DATA NeedNames/'GM1', 'GM2', 'GMB', 'GM4', 'GM5',&
        'GM6', 'GM7', 'GM8', 'GM9', 'EMRAT', 'GMS'/

    CHARACTER*6  AUname(1)
    DATA AUname /'AU'/
    !--------------------------------
    CALL selcon(NeedNames, 11, GM)
    GM(10) = GM(3) / (GM(10)+1.0d0) ! Moon
    GM(3) = GM(3) - GM(10)            ! Earth

    CALL selcon(AUname, 1, kmAE)
    !write (*,'(3(d30.18,1x))') GM, kmAE
end

